Intrinsic anisotropy of thermal conductance in graphene nanoribbons 
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Thermal conductance of graphene nanoribbons (GNRs) with the width varying from 0.5 to 35 nm is sys- 
tematically investigated using nonequilibrium Green's function method. Anisotropic thermal conductance is 
observed with the room temperature thermal conductance of zigzag GNRs up to ~30% larger than that of arm- 
chair GNRs. At room temperature, the anisotropy is found to disappear until the width is larger than 100 nm. 
This intrinsic anisotropy originate from different boundary condition at ribbon edges, and can be used to tune 
thermal conductance, which have important implications for the applications of GNRs in nanoelectronics and 
thermoelectricity. 



Graphene nanoribbon (GNR) is believed to be a promis- 
ing candidate for future nanoelectronic and spintronic de- 
vices i 1 ' 2 ' 3 ' 4 In comparison with intensive research on the elec- 
tronic transport, the study of the thermal transport in GNRs is 
not satisfactory until now, despite its importance for nanoelec- 
tronics. At present the chip-level power density in integrated 
circuit is on the order of 100 W/cm 2 , similar to that of nu- 
clear reactor. Increasing circuit density further would induce 
the exponential increase of power density, so thermal manage- 
ment on individual nanoscale devices becomes vital to ensure 
stable operation of integrated circuit. 5 This indicates a full un- 
derstanding of thermal transport property of GNRs is critical 
for developing any practical graphene-based devices. 

It was previously demonstrated that basic device building 
blocks can be constructed from GNRs with different widths 
and edge shapes, which dominates their electronic properties. 3 
For the evaluation of the overall performance of GNR devices, 
it is highly important to explore thermal transport characteris- 
tics of GNRs with different edge shapes and widths. In fact, 
a recent work found anisotropic thermal conductance in sil- 
icon nanowires (SiNWs): (110) SiNW exhibits a room tem- 
perature thermal conductance 50% and 75% larger than (100) 
and (111) SiNWs.— In contrast, systematical investigation on 
whether GNRs with different edge shapes give anisotropic 
thermal conductance is still lacking, despite recent progress in 
the calculations of thermal conductances of GNRs^ A clas- 
sical molecular dynamics simulation studied thermal conduc- 
tance for the 1.5 nm wide GNRs with different edges, 9 but 
the thermal conductance obtained obviously exceeds the up- 
per ballistic bounds.iSiil 

In this work, we systematically investigate thermal conduc- 
tance of GNRs with different widths and edge shapes, and find 
that thermal conductance of GNRs is primarily proportional 
to the ribbon width except for GNRs narrower than 2 nm, and 
more interestingly, exhibits strong anisotropy. Typically, ther- 
mal conductance of zigzag GNRs (ZGNRs) are 20% ~ 30% 
larger than that of armchair GNRs (AGNRs) at room temper- 
ature as the ribbon width ranges from 0.5 to 2 nm. Such in- 
trinsic anisotropy decreases with increasing width, but is still 
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FIG. 1 : Schematic illustration of two-probe transport system, where 
the center transport region (C) is connected to left (L) and right (R) 
semi-infinite thermal leads, (a) iV a -AGNR, (b) jV ; -ZGNR. 



larger than 10% even when the width is as large as 35 nm, 
and thus may have significant effect on the performance of the 
GNR devices. 

Nonequilibrium Green's function (NEGF) method is used 
to study thermal conductance of GNRs. This method is 
powerful to treat many body problems when interactions are 
weak, and can exactly deal with ballistic thermal transport, 
which gives the maximum thermal conductance of a mate- 
ria L 12 ' 13 ' 14 ' 15 Since the phonon mean free path is very long 
in graphene (~ 775 nm at room temperature) 16 , the thermal 
transport is nearly ballistic in pristine GNRs, and thus ballis- 
tic thermal conductance is studied in this work. First of all, the 
second-generation reactive empirical bond order (REBOII) 
potential, 17 which was proved to give phonon modes of GNRs 
compatible with density functional theory calculations,— is 
employed to calculate force constants. Then, Green's func- 
tion and the transmission function H(w) of phonon can be cal- 
culated for the two-probe transport system (shown in Fig.Q]), 
and finally thermal conductance cr is obtained by Landauer 
formula. 15 Since GNRs with larger width (W) will have more 
phonon transport channels and thus larger thermal conduc- 
tance, the scaled thermal conductance, defined as thermal con- 
ductance per unit area (cr/S ), is introduced to describe ther- 
mal transport properties for materials with different widths. 
Herein, the cross sectional area S is defined to be S = Wd, 
where 5 = 0.335 nm is chosen as the layer separation in 
graphite. Following the conventional notation, A^,-AGNR 
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FIG. 2: (Color online) (a) The scaled thermal conductance (cr/S) 
at 300 K versus width (W) for ZGNRs (black square), AGNRs (red 
circle), ACNTs (blue triangle) and ZCNTs (cyan diamond). The inset 
shows cr/S for ZGNRs and AGNRs with the width varying from 0.5 
to 35 nm. (b) The anisotropy factor {rf) versus width (IV) at 100 K 
(black square), 300 K (red circle) and 500 K (blue triangle). The 
lines are drawn to guide the eyes. 



(N ; -ZGNR) denotes an AGNR (a ZGNR) with N a (N : ) carbon 
dimer lines (zigzag carbon chains) across the ribbon width, as 
shown in Fig.Q] For the sake of comparison, armchair carbon 
nanotubes (ACNTs) and zigzag carbon nanotubes (ZCNTs) 
are also considered and their widths refer to their circumfer- 
ences. 

Fig. |2ja) shows the scaled thermal conductance (cr/S) at 
300 K versus width for GNRs. It is interesting to see that 
GNRs exhibit different properties with respect to carbon nan- 
otubes (CNTs). ACNTs have nearly the same scaled ther- 
mal conductance as ZCNTs, consistent with previous theo- 
retical results 11 . Moreover, the scaled thermal conductance 
of CNTs changes only slightly as the width changes. In one 
word, the scaled thermal conductance of CNTs is almost in- 
dependent of their chirality and size. This, however, is not the 
case in GNRs, where significant size effect and anisotropy of 
thermal conductance are observed. For GNRs narrower than 
2 nm, the scaled thermal conductance decreases rapidly as 
the width increases; while for wider GNRs, it changes very 



slowly. More important, ZGNRs have significantly larger 
scaled thermal conductance than AGNRs, indicating a strong 
anisotropic thermal transport in GNRs. For example, when 
the width is larger than 2 nm, the scaled thermal conductance 
of ZGNRs at 300 K is around 4.2 nW/K/nm 2 , which is close to 
that of CNTs, while that of AGNRs is about 3.6 nW/K/nm 2 . 
The room temperature thermal conductance per unit area of 
wide GNRs is of the same order as that of pristine SiNWs 
(~ 1 nW/K/nm 2 ) 6 . However, graphene has much larger room 
temperature thermal conductivity (~ 5 x 10 4 W/K/m)i2 than 
SiNWs (~ 6 W/K/m for a 22 nm diameter SiNW)2£ because 
the phonon mean free path is much longer in graphene than in 
SiNWs. 

To give a quantitative description of the anisotropic ther- 
mal conductance, we further define an anisotropy factor for 
the thermal conductance as r\ — [(0"/,S) Z gnr/(07S )agnr] ~~ 1- 
Fig. |2jb) shows the anisotropy factor rj of GNRs with the 
width ranging from 0.5 to 10 nm at different temperatures 
(100, 300 and 500 K). It can be seen that narrower GNRs gen- 
erally exhibit stronger anisotropy. However, the variation of 
the anisotropy factor is irregular at all temperatures when W 
is less than 2 nm. At 300 K, the anisotropy factor changes 
irregularly from 3 1 % to 20% when the width varies from 0.5 
nm to 2 nm, and decrease monotonously to 14% as the width 
increases to 10 nm. When the temperature changes from 100 
K to 500 K, the anisotropy factor changes little except for very 
narrow GNRs. Moreover, the anisotropy factor will decrease 
rapidly when the temperature decreases from 100 K (the data 
is not shown here). 

One can expect the anisotropy factor of GNRs will decrease 
to zero when the width is large enough, since thermal conduc- 
tance is isotropic in graphene sheets.— It is important to de- 
termine the critical width at which the anisotropy of thermal 
conductance disappears. We have extended our calculations 
to the ribbon width up to 35 nm. As shown in the inset of 
Fig. a), the scaled thermal conductance varies very slowly 
with the width for wide GNRs. The anisotropy factor of 35 
nm wide GNR still has a value of 13% at 300 K. Obviously, 
a direct calculation of thermal conductance at critical width 
using the NEGF method is beyond our computational capa- 
bility. So we apply linear regression to fit the data from 4 to 
35 nm for ZGNRs and AGNRs, and find that the anisotropy 
may disappear when W ~ 140 nm at 300 K. 

What is the origin of anisotropic thermal conductance in 
GNRs? It was once attributed to different phonon scatter- 
ing rates at the edges. 9 For ideal structures without defects, 
however, there will be no boundary scattering at all, since the 
boundary structure is incorporated into the phonon modes. In 
fact, our results show the anisotropy still exists even without 
edge scattering. It was previously revealed that the anisotropic 
thermal conductance of SiNWs can be reproduced from the 
anisotropic phonon structure of bulk silicon 6 . Differently, 
graphene, the bulk counterpart of GNRs, is isotropic in ther- 
mal conductance. This suggests that for GNRs, boundary con- 
dition at edges instead of the bulk property is the origin of 
anisotropy. 

The effect of the boundary condition on thermal conduc- 
tance is rather complicated. Since lattice vibration is a collec- 



tive behavior, almost every phonon modes will change if the 
boundary condition varies. Moreover, thermal conductance 
is always contributed by many phonon modes, and thus the 
total effect can not be demonstrated by focusing on partic- 
ular phonon modes, but should be analyzed from the whole 
transmission. Fig. Oa) shows the phonon transmission func- 
tion for 16-ZGNR and 28-AGNR, which have a similar width 
(W ~ 3.3 nm). Only negligible difference in transmission is 
observed at low frequency region (co < 100 cm -1 ), implying 
that the anisotropy disappears when the temperature is lower 
than 25 K. This indicates that the boundary condition has lit- 
tle influence on low frequency phonon modes. In contrast, 
obvious effect of boundary condition can be seen at high fre- 
quency region (especially at 400 ~ 600 cirT 1 and 1400 ~ 1650 
cm ): the transmission of 28-AGNR is always lower than 
that of 16-ZGNR. In principle, transmission function of bal- 
listic transport is determined by the number and the dispersion 
of phonon bands. More dispersive phonon bands would give 
larger transmission. Since the two GNRs have nearly the same 
number of atoms per unit volume, the discrepancy should not 
be caused by different number of phonon bands. Comparing 
with zigzag edge, armchair edge indeed gives more localized 
lattice vibrations and less dispersive phonon bands, resulting 
in the anisotropic thermal conductance in GNRs. 

The dependence of the scaled thermal conductance 
on temperature is illustrated in Fig. [3jb) for 16-ZGNR 
and 28-AGNR. As expected, the conductance increases 
monotonously with increasing temperature. Evidently, 16- 
ZGNR always has larger scaled thermal conductance than 28- 
AGNR. The anisotropy factor is small at low temperatures and 
keeps around 16% from 100 K to 500 K (data is not shown). 
Since the scaled thermal conductance only slightly varies with 
the width when the width is larger than 2 nm, it can be used to 
determine the maximum phonon thermal conductance at dif- 
ferent temperatures for wide GNRs. Thermal conductance of 
a ZGNR (an AGNR) (W > 2 nm) can be estimated through 
multiplying the scaled thermal conductance of 16-ZGNR (28- 
AGNR) by the cross sectional area. 

To further understand the effect of boundary condition on 
anisotropic thermal conductance, we artificially impose fixed 
boundary condition on the above two GNRs (i.e., 16-ZGNR 
and 28-AGNR) by setting the mass of edge atoms to 1.0 X 10 7 
atomic mass unit in the calculation. Fixing edge atoms causes 
a decrease of scaled thermal conductance for both GNRs, as 
shown in Fig. Ob). The scaled thermal conductance of 16- 
ZGNR decreases noticeably, especially at high temperatures, 
whereas the decrease of 28-AGNR is relatively small. Af- 
ter imposing fixed boundary condition the anisotropy becomes 
smaller but does not disappear. This shows the anisotropy is 
robust even in this extreme condition, which is substantially 
different from the free or periodic boundary condition gener- 
ally used. 

It should be noted that both phonon and electron transport 
can contribute to thermal conductance in practice. The elec- 
tronic thermal conductance of gated GNRs was well discussed 
in Ref. |22j. Since recent experiments showed that all sub- 
10-nm GNRs are semiconducting, 23 thermal conductance in 
GNRs narrower than 10 nm should be mostly dominated by 
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FIG. 3: (Color online) (a) Transmission function versus phonon fre- 
quency (a>) and (b) the scaled thermal conductance (cr/S ) versus tem- 
perature (T) for 16-ZGNR (black solid line) and 28-AGNR (red dot- 
ted line). The inset shows the decrease of scaled thermal conductance 
(A) induced by fixing edge atoms. 



phonons at all temperatures. While for GNRs of several tens 
of nm wide, the contribution from electrons might become 
important at low temperatures. 

In conclusion, we have systematically studied thermal con- 
ductance for ZGNRs and AGNRs with the ribbon width rang- 
ing from 0.5 to 35 nm. Similar to electronic property, thermal 
conductance can also be tuned by the ribbon width and edge 
shape. The scaled thermal conductance {cr/S) decreases as the 
width increases and becomes nearly width independent when 
width is larger than 2 nm. More interestingly, ZGNRs have 
13% ~ 31% larger thermal conductance than AGNRs at 300 
K, showing that thermal conductance is anisotropic in GNRs. 
Unlike in SiNWs, the anisotropy in GNRs is caused by dif- 
ferent boundary condition at edge but not the bulk phonon 
structure. Comparing with zigzag edge, armchair edge gives 
flatter phonon bands and behaves more like a fixed boundary. 
The intrinsic anisotropy found here is not limited to GNRs, 
but is also expected to appear in other materials with hexag- 



onal structure, and is strong enough to be observable experi- 
mentally. Our work provides useful insight into the thermal 
management in GNR-based nanoelectronic and thermoelec- 
tric devices. 
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